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Zqua'^ions of transfer v/hich include, the effects of therr.; ! enission 
end anisotropic muitipie scenterln^ are solved ^ by the method of discrete 
ordinates., and the angular distc„.^u cion- and net. flux of ..outgoing radiation 
are nlculated for various model, clouds and clohdy atmospheres . The 
elo'.’d models addp.ted ara piaiie-paraiiel and of normal optica.l thickness 
T" 52 0,1, 0.6, and co, . The principal wavenuniber. intervals chosan for 
s-xamination are-;cdncare^ about 889,5 cm“l and 1173.5 cm"!. Th%.corstrii» 

■ bu cions of the. outgoing radiation field at -the- top of .. the .atmosphere, , 

'due.ito diffuse reflection,' diffuse transmission,, direct transmission, dnd 
' thermal . emission are traced -out in detail for. each model clbu.dy atmosphere. 
Based on the expllcic outgoing raiiation field predicted for each rao'del, .• i 
the following four; general conclusions 'are reached: -(1> limb darkening s| 
is- considerable at very large zenith angles, and. especial ly.prpnounced . I 
for- thin clouds -which ace.nych cooler than the effective redialing j 

iayfer bdlcw them, (2")- A considerable' flux surplus’ in= the iiiffared ' ,j 
window compared with the 'rest of- .the infrared spectrum is predicted ^fbr ■ 
atmospheres- containing cold, ^.thin cloud’s. <3)^ 'A jabdest flux dericife in / 
■the infrared wihdox'j is predictad for' atmospheres containing very .extensive/ ' 
warm, thick -clot ds. (4> For the same cloud- top temperature, the. butfoingj 
-radiation fields are virtually t'ne' san® for an isothermal cloud and a ' j- 
cloud-having a realistic temperature jgradient^ of 5®C km“l. Ail -con- . 
elusions have to be modified if iaultiy^i,^loud decks and cioud eonfigu- 
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1. Introduction 


A knowledge of the transfer of infrared radiation in and 
through clouds is of fundamental importance in the overall 
investigation of radiative transfer in the Earth’s atmosphere . 

It will '>e the purpose of this paper to further progr-ess toward 
placing the physical understanding of the radiative character- 
istics of water clouds on a more nearly equal footing with that 
of the clear atmosphere. To this end we shall be concerned 
with: a) the numerical determination of the role and importance 

of diffuse reflection and transmission and of thermal emission 
in and from clouds, b) the angular distribution and net flux 
of outgoing radiation in selected wavenumber intervals at the 
top of certain model cloudy atmospheres, c) an interpretation 
of the main radiative features in terms of the relevant physical 
parameters, and d) a qualitative extrapolation to probable 
radiative characteristics of clouds and cloudy atmospheres at 
other regions of the spectrum and for other physical situations 
not explicitly covered by the computations in this paper. 

In order $o maintain physical rigor it is imperative to 
eliminate approximations based on intuition insofar as possible 
in the mathematical fomulation of the general problem. It is 
further required to maintain as much generality in the formu- 
lation as possible in order that the range of validity not be 
overly restrictive. These objectives require solutions to 
rather sophisticated equations of radiative transfer; a powerful 
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Lochaique developed primarily by Chandra, -hhar for tlic purpos-.c 
oj obtaining these solutions will be the one exploited here, 

A full account of the Mie theory for single scattering and the 
theory of the transfv^r of radiation through the surroundin-g 
ga.seous atmosphere will be included in the computations which 
are .required in order to fix the boundary conditions and des- 
cribe the microscopic scattering processes. A complete des- 
cription of multiple scattering processes then follows from 
the solutions to the “relevant equations of radiative transfer - 

2. The Equation of Transfer 

The cloud models adopted are plane-parallel and of arbitrary 
normal optical thickness. The particle size distribution is 
maintained Independent of optical depth, although the mass dis- . 
tribution with height is left arbitrary. Each mass element in 
the cloud is assumed to emit radiation thermally in accordance 
with Kirchhoff's law and scatter radiation in accordance with 
the Mie theory. The temperature is assumed to be a monotonic 
function of the optical depth and Independent of horizontal 
displacement. The (thermal) radiation field is assumed to have 
axial symmetry about any normal to the plane of stratification; 
this immediately eliminates any outside discrete source of 
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I’r.ciiation (e.g. the sun) from consideration. The state of 
polarization of the radiation field is neglected, thereby 
roe cing the transfer equation to a scalar equation. Since 
v/o are interested in both the magnitude and angular distri- 
; bution Ox radiation, the specific intensity is our choice 

I 

I 

I for the dependent variable. In order to maintain as much phy- 

sical rigor as possible, no quantities which are not strictly 
independent of the radiation field are allowed to enter the 
equation as free parameters. 

It has been shown by Samuelson (1964) that the basic eqi tion 
of transfer which satisfies these requirements is of the form 



- (/- ia) ^ 


( 2 . 1 ) 


where the relevant parameters are defined below. 
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The specific intensity is the rate at which 

radiant energy confined to a unit solid angle and a unit wave- 
nu.r.’jor interval crosses (at a level f ) a unit surface area 
normal to the direction of propagation. The direction cosine /U 
specifies the cosine of the (zenith) angle between the direction 
of propagation and the outward normal to the plane of strati- 
fication. 

The normal optical depth T is defined by 

T =■ ^ J 
2 . 

Where is the heii^t of ttie cloud ^op, a is the height of the 
level of interest and is the total number of particles 
per unit volume. The effective cross-sections per particle for 
extinction ( Xg) , absorption ( ^^ > , and scattering ( ^/) asre 
defined by ' 





( 2 . 3 ) 


where ^ are respectively the Mie efficiency 

factors for extinction, absorption, and scattering for particles 
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of radius r. D(r) is the ratio of the number of particles per 
unit volume per unit radius range (centered about r) to A4 , 
i.o. the part’ le size distribution function. 

The Integrated phase function for single scattering 
is defined by 







[Qs nr 

g 

J QsVf'Di") 


\d¥'- f} ^ 


(2.4) 


wliere ^ is phase fuaction for single scattering 

througbi the angle d for particles of radius r« It can be sl^wn 
fr(^ spherical trii^moMtry ttot 41 nay be expressed in terms 
of the coordinate angles #, $*, 4 ^ throui^ the relation 


tdO- 0 ^ ^ (2,5) 


where s ce». and ^ r tea ^ are respectively the 
direction odsinos of the incident and scattered radiation referred 
to the outward nonuil to the plmie of stratifioati<m» and ^ and 
(if are tvhe oorTespoading aximnthal angles measured in the 
plane of strati fioatiom* 
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is normalized to the albedo for single 
scattering, , defined by 



J rx^) Jr' 

f dlgirr' £)(n 
0 


( 2 . 6 ) 


The quantities , and Cc^ 0} may be calculated 

from the Mie theory. 

Bi V" ) f the Planck function in intensity units, is assumed 
to be a monotonic function of T 

: It sfaO|ild.be noted that all "quantities in (2«i) are 

functions Of the wavehumher ; thus. (2,1) is valid only over a 
wavenumber interval m>3Xt, enough such .that the derived quantities 
remain sensibly constaht across the intex^al. It should further 
he noted at this poia' ':hat JCt^/4) refers the total 
radiation field, and thus implicitly includes all radiation frcmi 
the ground and bounding atmosphere which h^rreeii transmitted 
directly to the level T without any intermediate 

scattering and/or absorption and emissim processes. 
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3 Input Data 

Certain input data are required before the paran» 3 ter s , 

, and can be calculated. These data Include 

the complex index of refraction of water, , the 

particle size distribution, D(r) , and the wavenumber i'itervals 
of interest, , 

Figure 1 is a plot of the real i n) and imaginary (>c ) parts 
of n vs, wavenumber V • The curves are fitted to the data 
publiehed by Centano, •Johnson and Terrell, Herman, and Stephens 
as listed 1 :^ Pehndorf CI963) • The section of the A vs. -i/ " 
for the range (1050 cuT^^ 1/ 1450 was obtained I 9 ' 

rather arbitrarily weighting the more recent data by Herman and 
Stephens twice as he|^^ as the older data* This is a very 
critical section of curve, since only a modest shift of the 
value of H, will affect considerably the relative amount of 
absorption calculated'* 

TlM shaded iK^ticUis /^ tl^ fi^e i^^ 
intervals which will be considered in this paper* Table 1 gives 
the ran^d (wiNn«^ ^ ^ and 1 ^ mean JC^ * ) of these. 


in'^ervals in inverse centimeters* 



Table 1. 

V/avenumber /ntervals 

( ^ V < 

) and 


wean (^ ) used in the 

computations 

IT 


V-mivu 


V 

Interval 

(cm*"l) 

(cm“l) 

(cm~l) 

1 

475 

500 

487.5 

2 - 

871 

908 

889.5 

3 

1149 

1198 

1173.5 

4 

1425 

1450 

1437.5 


The first two :4ntervals were selected to be representative of the 
highest and .lowest yalues of n. The last two values were 
selected \f or the'. p testing the region of the spectrum 

w^ere % 'is a minimum i 


Tliie particle size distributions considered are taken (with 


minpr':ii;odifications) from Nei.burger (1949) and are represented 
in figures 2a-2d, The narrower distribution, p-j^(r), centered 
hbout , d, % is representative of California' stratus, 

excluding'' the bj^se of the cloud; the broader distribution, , 02 ( 1 ’) 
is representative of the base , 

tj. - , . ' ' ' ' ... 


fhe graphs of 


(2e and . ate superposed over D(r) . 


These, values., as well as all 'subsequent Mle scattering para- 
meters, were obtained from a progra^i furnished to the author by. 
i« i^nn and t* Michels of . the Theoretical Division, Goddard 
Spkcd I^ilght Centef* The actual computations were performed on 
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the IBM 7094 electronic computer. All minor wi^£fles in the 
calculated curves have been suppressed. 

Four points should be mentioned: a) Decreasing the wavenumber 

has the effect of spreading the curves out to the right, 
b) increasing tc has the effect of increasing , c) increasing 
% has the effect of increasing , and d) decreasing k 

beyond a certain point has the effect of magnifying the oscillatox'y 
nature of ~ • 

In general every mass element in a cloud will extinguish 
(absorb and/or scatter) incident radiation at a rate in accord- 
ance with an effective cross-section for extinction which depends 


upon the absorption properties of the pervading gaseous atmos- 
phere as well as the absorption and scattering properties of 
the particles present. Define g to be the probability that the 
incident radiation which is extinguished is extinguished by the 


pai*ticles« It can be shown (Samuelson, 1964) that, if g is 
i’^dependent of t , equations (2,1) - (2.2) remain valid simply 


by replacing tDa with ^ and Ai with ^ Ap 

Deviations of g from unity in the wavenumber intervals 
centered about V “ 487.5 cm“^ and ^ « 1437.5 cm*"^ are due 


almost exclusively to, absorption by water vapor. The exact 
amounts of these deviations depend very critically upon the 
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absolute humidity. To further complicate matters the water 
vapor spectrum is extremely variable over very small intervals 
of wavenumber, and the transmission functions over these 
intervals are correspondingly far removed from simple expon- 
entials. Hence, therefore also 

functions of 7" , and the equation of transfer cannot be treated 
correctly by the method used in this paper. 

In order to obtain an approximate idea of how water vapor 
affects , transmission functions (, T ) were calculated for 
various thicknesses of atmosphere containing saturation mixing 
ratios of water vapor. The two sets of values of the mean 
pressure and temperature used were, [P = 950 mb; T =» 287 Kj and 
[P = 417 mb; T = 259 K] and the increments of thickness of all 
layers were restricted to those which correspond to pressure 
increments of between one and five millibars. , An attempt was 
then made to determine an empirical absorption coefficient 
for each wavenumber interval considered, such that the expression 
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i 


is approximately valid for short path lengths and for the two 
wavenumber intervals of interest. The coefficient 
was found to vary by a factor of about three over the path 
lengths considered. Since the generalized absorption coefficient , 
( ) of Elsasser for water vapor, as adopted by Wark ^ al . 

(1962) for the same two wavenumber intervals considered here, 
ax’e v/ithin CD. the relevant ranges of , it was decided to 

replace with these values of explicitly in order to 

resolve any «'mbiguity 'in the choices of . 

Table 2 gives the values of X£ t ^ each 

of the particle . size distributions depicted in Figure 2. The 
inherent crudeness in calculating the values of ^ , requires 
that these values be considered only for the purpose of illus- 
tration. We shall therefore restrict ourselves for the raost|v- 
part lo the wavenumber intervals centered about v = 889-2^::^m 
and ^ = 1173.5 cm“^, where water vapor plays a negligible' 

i - 

role in the transfer problem. It should be noted that clouds 
will be most transparent in the infrared window, since elsewhere 
water vapor, as well as carbon dioxide and perhaps ozone, will 
tend to increase the effective, normal optical thickness 
[eq. (2.2); cf. also Table 2 and Figure 1]. 

f - - ^ - 

I - 

i - ■ ' ^ , 



Table 2, Tni:egrate<J scattering and absorption parameters characteristic of the 
particle ^ize ^{istributions and H 2 O inixing ratio considered. 
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The solid curves in Figures 3a-3d are the phase functions 
for single scattering integrated over all particle sizes and 
nor-nalized to ^ ^ » Thus 


J Qs Oir) >py.(s^ 

p(ura (9) z: 

J d'r 

and 


(3,1) 



(3,2) 


^yhere the integration in (3,2) is performed over all solid angles. 
The dashed curves in the figure are 14— terra Legendre polynominal 
expansions v/hich represent be&t fits in the sense of least 
squares. The fits are expressed in the form 


^ (3,3) 

where the coefficients are constants independent of & . 

These approximations to the relevant phase functions will be the 
ones used in the remaining numerical solutions. Although the 
fit for y = 487,5 cm"“^ is very good, not enough terms are 
retained to give completely satisfactory results at lu-•g'^r wave- 
numbers. Horo will ba said about limitations this will impose 
on the accuracy of the results in Section 8, 
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4 . Formal Solutions 

The specific intensity refers to the intensity 

of radiation in the direction yu at a. level f , In order to 
specify correctly the boundary conditions in a simple manner we 
shall require that refer only to the diffuse radiation 

field which has arisen • through one or more scattering and/or 
emission processes within the cloud itself ; thus any radiation 
from an outside source directly transmitted to the level Z will 
not be included in the- solution. Since (2.1) refers to the total 
radiation field, it must be modified accordingly. 

The intensity of the diffuse radiation field in the 
direction v/4 at a level z may be thought to be composed of 
three components: 

1) That component of the intensity which has arisen as 

a result of the radiation field (from the upper 
hemisphere) being diffusely reflected by the cloud 
into the direction . 

2) That component of the intensity which has arisen as 
a result of the radidtion field (from the lower 
hemisphere) being diffusely transmitted by the 
cloud into the direction yu , 

3) That component of the intensity which has arisen as 
a result of thermal emission by the particles within 
the cloud itself. 
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w'e sl:ci.ll suppose that the individual radiation fields do no 
interfere, and hence may be treated independently , 

Consider now a point source removed to infinity in the 

- 

out'..'ard hemisphere. Let be ‘che flux of beam radiation 

fi-cm this point source crossing a unit sux’face area normal to 
the beam; further let (r , ‘'y.) be the direction of propagation 
of this beam. The azimuth -dependent equation of transfer 
describing the resultant diffuse radiation field may be written 
as (ci. Chandrasolchar, 1960; Samuelson, 1964) 



*\Ie cons^ider the atmosphere outside the cloud, and the ground, 
to conotribute a radiation field which is made up of a discrete 
number^f point sources. The separa.te conti’ibutions from each 
point source are added up later in calculating the total 
contribution. 
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wiicv' phpvSe :ii;nc';:ion for radintion singly 

licaiior'^d ibro.igh the ar^gic defined by the direction of incidence 


/ V ' >' 


and the direction of scattering (/U 6) „ For a 


f urn eiori e:.'p:/essxole as a finx'co senes expansxon ox 
l.cg-endre poiync-r.ials [cf. eq. (3.3) j the solution to (4.1) can be 
forinally ’.vritten as 






(4.2) 


In particulax', the azimuth-independent term in (4.2) obeys the 


equation 


i 7^^/' 'A 

'A J, Ct-./^'O 




- i f 
-/ 


(4.3) 



r*zi;uuv.h-indcpendei;t toinus in 




''{/■.jy.) 'CO the 


4 . 1 i \ 


c: pans ion in sphooical harucr.ics of p, Cc^-s. 3>J in (3,3); 










Ay \7C so.all preset, tly see, the solutions for the renaining 
tor..:s in (d.2) are of no interest in the context of an axially- 
syia-etric radiation f:'reld, v/e should further note at this time 
jhat (2,1) is the correct equation governing the diffuse radiation 
field arising frori: thex-mal emission; i.e., no laodif ication of 
(2,1) is I'equired in describing the third intensity component 


cited at the beginning of this section. 

The method of discrete ordinates (Chandrasekhar, 1930; 
Samuelson, 1934) is the method used in this paper, to solve (2,1) 
and (4.3). The continuous radiation field is replaced by 2n 
linearly independent beams of radiation^ n each in the upward and 
downward directions. This artifice allows the integrals in (2,1) 
and (4.3) to be replaced by sums, and the resultant system of 
linear first order differential equations admit closed solutions. 
The solutions are not analytic in , of course, since they ai»e 
evaluated only at 2u discrete intervals. 
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x'he solution to (4,3), in the ntli ap:)i-o::imition , is of 
tho form 


’crA-; ^ jFj Z 


^ Atcs'"'' '’"V -t - / ^ \P^ 1 




-Z Jj 

i_>2 = 0 


,_J 


n Z 

rf 

L T , 

.cT; 


p y 

+ i 

i r 

li-A-:^,, i 





(4.5) 


± ±n) 


) 


where 


= - 


2£^l- 3, 


t/ 

1 

•*o 


A (JL*0 ^ 


? 


— - f 1 

j 


0 

I 


y 


(ji- ^ A/-0 




(4.6) 


4/ 


'=!?“!.■ 




J>s±, 


“X 

7f/«A 


(/i ±/^...^±»)^ Cf.7) 
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H&) - 


jf 

X = t -- 


(4*8) 


' ' 


II 


Once n and K .*ave been chosen, the discrete, in ter?a*»s i 

{i . - : , * ' i; 

,,,, ±?x) and T:he Gaussian weights A/' Ca^ 

o' ‘ V. , 

may be obtain- 1 respectively as,, ,the zeros of ,the Legendre poly«- 
nomial RyJyM ) formula ? . 




1 


4u j ), ^ : 


(4.e) 


it is always true that 


A* - V'a' 


-a,- = A/ 


u.io: 


, -- 
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foi* all values of n. For X = 13, n " 7 ( che values used in 
this paper) the values of ;:ind ‘ ' j 

found to be those listed belo-v in Table 3. Equations (4.6) and 
(4.7) can then be solved interdcpendently on the electronic 
computer to obtain the various values of f.Ot) 

This leaves only the 2:i constants of integration M;te< I 

to be solved. The requirement that the diffuse radiation field 
contribute nothing from the upper and iov/er hemisphere at 
respectively the top and the bottom of the cloud yields the.2n 
boundary conditions 






(4.11) 


From (4.11) we obtain the 2n equations required to solve for the 
2n unknowns, M±a. (oi ' «>: (1^. The angular dependence of 

the. outgoing intensity at the top of the cloud is then given by 
' , , and, at the bottom of the cloud by 

(i=l, 7). If the cloud is semi-infinite in extent, all the 

terms in (4,5) containing > -j are suppressed, 

and only the n conditions ' 

I - 6 (Zz 7^ (4.12) 

ajje applied. 
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'•.••c.bls 3 


(- 


ssin,n 3i^:hts ana ..i 
l-'I-Poixi'C .’uadratui-e 
-:-l) 


Crete .'ntervals 
vor;nulii over the 


eievant 

/nterval 


# 

f 


dP - 

.1 

A O Q C O Q O 1 

0 o 33511946 

2 

A O O 

V> o O*-- O 

0.03015809 

2 

0.S2720132 

0.12151S57 

4 

0.6S72&290 

0.15720317 

5 

0.5'l52<SS4 

0.18553840 

6 

0,31911237 

0.20519S46 

rt 

i 

0,10805495 

0,21528385 


Let be the axially-symmetric intensity of 

radiation v/hich has arisen as thermal emission in the outward 

j , - 

hemisphere and is incident on the cloud at a level T and in 
the direction (formerly )o It is convenient to 

define a scattering function S( ') and a transmission 

function T( Z; j j ^ such that the./ intensity of radiation 
from the entire upward hemisphere which is diffusely reflected 
by the cloud into the direction ^ is given by 

i 

^ [sCt, 


( 4 . 13 ) 
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:.rid tii3 intensity ot radiation 2x'oiu the ortire lower hemisphere 
diffusely transmitted through the cloud into the direction is 

given by 


/ 



(4„14) 


where and 2~=^ are respectively the normal optical depths 

at which the cloud bottom and the .cloud top are located. In the 
context of the solution to (4.1) it can be shown (Sarauelson, 1964) 
that S and T obey the relations 


^ <4.15) 


and 





2rr 

f 

O 


(4.16) 


or, by virtue of the axial synanetry eiihibited by and the 
form of 





(4.17) 
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ind 








'■ V^T 




(4.18) 




rl-.-M'! 


uld be ncl.j_ biir.t tbe sense ox direction of CT/ is 

reversed in the present conte;-:t from what i'C is normally. 

Srnations (4„13) and (4.14) comprise the solutions to the 
components of the total outgoing intensity at the cloud top v/hich 
arise in the first case ^ron the upper, and -in the secona case 
from the lower, hemisphere, and are respectively diffusely 
reflected from and transmitted through the cloud. Upon, replacing 
B( T ) with the finite power series representation 




B&-) - 


IT 


(4.19) 




it can be shown (Samuelson, 19S4) that the solution to (2.1) for 
the radiation thermally emitted by the cloud is of the foimi 


r. 










aJ 


^ r-fj 

, r-0 ^ 


y;s 


r? 





J 


( 4 . 20 ) 
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I’o tho values of can be generated by neans of the 




r,s 

:us recursion relations: 


a, y 

.■yJ 








» 






// 


r 

^rjM 


i C:/-r) \ 

-/J 

i _ I 

‘ XCf'^}-r\ J 

UrfO^rjj) 

L;2(://'-r> - 3 J 


I 


y>/\ 


A/ 








'J 


■r, r 


— lO^^ T- 


rt 

) ' 


^ ” O'; hrJl 


and 


(r= Oy..j M-x) 


) 


(4.21) 

(4.22) 

(4.23) 

(4.24) 

(4.25) 
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* I t 




) - 





S-'r 



(4.26) 




i t •. t * i j ^ 


3 ’--j "'^'■'^3 . 


'21i 3 oroccdurc lov solving (<I-o20) is exactly the saaue as that ior 
(<!.. 5); in particular the boundary conditions (4.11) are incor- 
porated in solving ior_ the 2n cons C C'O'i-Ji C<^~ /j •‘ •j » 

it the cloud is semi-infinite in extent, conditions (4,12) are 
used, and the constants M~d C<- •■•j 

all supnressed. 


5. Cloud Top Radiances 

C "sider novi a plane-parallel model cloudy atmosphere composed 
of an optically finite cloud imbedded in an otherv/ise non- • 
scattering atmosphere, bounded on the bottom side by the ground, 
and on the top side by outer space. Let ( ^ t' 2 ^ ) 

be respectively the optical depths at v/hich the cloud top , cloud 
bottom, and the ground are located. Further let 



( 5 . 1 ) 


i 
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L w Ox V i - -j X -* L- j . i, t- y ^ 

rhe lov/er 

vliich is dii-^ctly ii’ 

ansnitted 

,r._ 0.-1 (..'JilS- , >;j 

froin the 

yl'' f is 'ciiOu 'uaIC 3 Uiii Oil 

the coinpon 

:rans:::ittoi thi*cu;;h, aiiu 

diffusely 


tlioi’./.ally enittod from the clcucl„ Thus 




J 


(3.2) 


or, by virtue of (4.13), (4.14), (4.20), and (5.1). 
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•^Vi'.iuatioa ox the in v 0 2 ,'xCvl,-i in (o,3) requires a 
su.v.'.xation . Xn order to hold the niiv.iber ot terms 

.vhiio still main /nxnin;-; accuracy, the Sykes double- 
dr'.tero xormula ('-inj; and -'loronce, 19S4) is the one 
vofjrals are appro;:iv:'.t-^d by 
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o Ca^U 
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'■'■a) jV-'';) — . \/^‘) ) ^yy‘ , (Sol) 
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v'here the v/ei.qhts and the discrete intervals yu,- are listed 


i 


in Table 4 . 


;le 4« -deights and discrete intervals relevant to a 
4-, dint r aadrature rci-mula. xver the interval 

(0, -rl) " 


2 

3 




0„93055316 
0.66993052 
0. 3300094 S 
0.06343184 


0.347854G 
0.6521452 
0 9 65214 o2 
0.34 7854 8 





'■..'lie ...i.;,.osiihere chosen vo voprc;: e.'.v. v.>^ilcr.i concil'cions 
•j CnliJ.'o-.-uia stratus r.ii{i'ht los;i.c:.lly occur is a'urivosphore ir55 
T:;' '..'avl:, Yar.iar.io to , and Llonescli (l£Gi) and reproduced 
.ii:-r r.rCdiiications in Taolo d i,eio-/„ Ihe intensity components 
£, /• ' and Y-; . •-* = v/cro cLtained from 

the prassure levels ? = 417 mb and Y = S5C mb frora an IL.l 7094 
computer prograia for calculatin,-]; the infrared radiance of a 
spher.'i.cal atmosphere which was furnished to the autlioi- by 
IP. Van Cleef of the ileteoroloyical Satellite Laboratory, U. S. 
V/eather Bureau (ci. V/ark, Alishouse, and Yamamoto, 1SS4) . 
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'"he c:l ;j.d r.iodolo studied •. xll ^ostricted to tln.-oc 
.,1 , :.iv;kaos3Cs ( 4 = Col, C.5. c- \.\.o cloud 'cop p-"cssurs 
u.- , :_Totuiu; levels (? - dl.? n'r/. 'd = 3c9 I.! sad - 950 mb, 

-- „17 uad one particle s_se distribution „ It v/ill 

So.; a in Section 9 tnat realistic c-vioUtions from isothermal 
cloud .models produce noi;li^ibl:' chances in the computed outgoing ' 
Intensities o Fir.ing the optical thickness oi isothermal clouds 
rencvjrs a description of the angular distribution and net flux 
of outgoing radiation independent of the actual mass distribution 
of rater in the models [cf. eq. C2«2)]o • Ihe two wavenumoer 
intervals considered in detail are the ones (Table 1) centered 
about V = 889.5 cm”"^ and v == 1173.5 cm“l. 

Figures show the cemppted i’esuits for the wavenumber 

interval centered about i/ ~ SIS, 5 cm*"^. The dotted curve (5) 
in each separate figure refers to the diffusely reflected com- 
ponent of the intensity [cf. eq. (5,3) 


/ 


while the lov/er solid curve (E) refers to the thermally emitted 
compoaeat, 


I. . 

^ y 



n:-.. .'-.l?uy accouutod i'or j ^,/u co;i;putr.'.:ions cf the 
cc,„ L'l'.e t.;o curves rofei" to the terus 


:. O'' 




i i ’i''' ^ ; /"v ' “ ■■ 1- ' ^ ' ;- 0 






iudicu'.ed in the figures by (D) and (T) respec tively, the 
ccivi," oaentL ox the intensity v/hich originate in the loyer he;ais-' 
phsre and are ’espectively directly transmitted and' dlfiiiseiy 


trrinsnitted through the cloud into the direction yOL , These 


last t’.vo components, of course , do not contribute to the net 
radiation field if the cloud is semr-infinlte in extent. 

The upper solid curve in each figure is the sum CO of the 
four corresponding components, and thus reefers to the net 


radiation fid 


U tt 


The dash-dot curve in each case is the standard 


of comparison (3', i.e., the ter al intensity v/hich would result 
at the cloud top if tiie cloud were of the same optical thiCtotes;S|/ 
but composed of nonscattering particles, computed from the 










expression 




-SCs)]i£^ , (5.53i 
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v/hoi'j BCI'', ) is ihe r-lruick xunv.'i;ion i.i in j^nsity units integrated 
over tjie r; levant v/avenumber interval and co:nputed for a cloud 

xrri'oorat'are T „ 
c 

"nrhans the racst reriiarhacie general feature evident in 
r'ijiercs 4n-4i is the lack of -,a overall narked de-viation fr-om 
greynkss fqr the realistic clcuci models^ in spite ox the fact 


tlae albedo fOx single scrtterinr is ra'i 


This is especially ti’ue of xhe lev/ cloud models, v/here the 
difference in teraperatjare bat^veen the cloud and ground is 
- . ' ‘ relatively small » 

i Figures 5a-5f are the same as Figures 4a-... except that the 

4 - •• - 

4' T V y/aveaumber interval is nov/ centered about *P = 1173.5 cm’'^, and 

^ Is the higher value of .76S. The differences betv/een the 

4-:-=ns.-- ' 1 

|1 1 upward intensities relevant to the realistic cloud models and 

< those of the corresponding standards of comparison are more 
f marked than the differences illustrated in Figures 4a-4f. This 

4 ■ _ is ^clearly a result of the higher albedo for single scattering 

. ,in uc v/avenumber interval centered around t 7 = 1173 ..5 cm”^. 

f ' \ In both \7aveuuinber intervals there is a surplus of radiation 

with the standards) for optically thin clouds at 
t b - ^smail zenith angles and a tendency toward a deficit of radiation 


If . 


% '''■ rl- 


i;k’ ' 







o o 

— o — 


.nj;los« U_:on i.:‘; inv -vco.upariscn O" ?i--ures 4a-4d 


.'r.-od iOv Si.iai! 


j iv. is noticed that the 


r. L ind.: ci the serplns inei-oasos as Cf . la.-^ ^ and the 


v.torencs (T^-T ) tet^.•eon -che 
^ o 


rrc-und and the cloud 


. 'je- asc; \,„e i'nnctional dependence seen /-e, , ben'ever, stays 
.tc./: nhc; sa:..e, Apain, upon interco;uparing Figures de- 4 f and 
^ -2»'- i« ci^at «^nis surpass CwC cj^ie 3 a ao x x ci o as 


and the magnitude oj 


T ~r ,** - C! /* o 


licit increases as 


increases. It is quite clear, therefore , that the sui'plus 
of radiation is primarily due to the radiation diffusely trans- 
mitted through the clouds; this radiation is. lost from the 
radiation field in the cases of the ' standards- of comparison, 

V.'e must look elsewhere for explanations for the form of 
the angular distribution of radiation and the deficit of . radiation 
from thick clouds. These explanations are readily obtained 
from physical considerations discussed in the next section. 


C, Pnysical Considerations 

Figures 6 and 7 illustrate the effe ■ ’.ye emissivities of 
semi-^nfini te model clouds for the wavenumber intervals centered 
respectively about y = 889,5 cm“^. and v = 1173.5 cm”^ as a 
function of , In either figure the upper- solid curve (AS) is 



the c.. ‘tec Vive emissivity L€ ~ i£,CT-)^ ox a cloud having 

a 2>av.ic*-j size distribution 0,C"J ^ Tiie lov/er solid 

cvrvo '’IS) in each iigure illustrates the eifectivs enissivity 
oi a clovl ccr.v'osed of particles scattering radiation isotropically 
the sar.io respective values of . These latter curves 

were obtained v/ith the aid of the expressioii (cf. tin. Appendix) 




> 


( 6 . 1 ) 


v/here A~i/ is the wavenumber interval of interest and H(/x) is 
Chandrasekharas H-function (cf. Chandrasekhar , ISSO; Table XX, 
Chapter V) defined by the nonlinear integral- equation 


/^i/) - / z^o/^ f -^^0 . (6.2) 


o 


It is immediately evident from the figures that the fox'ward 
scattering nature of realistic particles is responsible for .the 
higlisi’ values of g , and that higher values , of 60^, tend to 
reduce the -magnitude- of £ in all cases. 

- Consider now a semi-infinite’ plane-parallel isothermal 
cloud composed of partly thermally emitting, partly isotropically 
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scci'cti; rir.y _:,',v'i.icio3* -■-j.l ''cl:u oarviclus cniit radiation at the 

3ar;.i' rate they arc all ac the same tenperarure „ Kov/ever, 

the pa:-ticlos at the top of th j cloud v/iii scatter I’adiation 
at a rate of only about half, oay, of oo.o rate at v/hich particles 
d:-p v.-ithin the cloud will scc/cter radiation, since only about 
half of the radiation available to the particles deeper in the 
cloud is avr.ilrble to the particles near the top (the upward 
hemisphere cu.n contribute nothing at the top) , Since the 
radiation is assumed to be scattered isotropically by each 
pai’ticle,- the rate at v/hich radiation is scattered by any particle 
into any direction in dependent only upon the level T where' the 
scattering takes place. Hence, the intensity of outgoing f 

radiation at the top of the cloud will be enhanced in the direction 
- 1 over the intensity of radiation in the direction O , 

bi-cause deep laj^ers contribute more to the outgoing intensity in 
the former case than in the latter for equivalent optical path- 
lengths. Limb darkening is thus explained on physical grounds. 

We now overlay the cloud, v/xth an optically semi-infinite . 
isothermal slab held at the same temperature as the cloud and 
allow the radiation field bet-veen the two boundaries to 
achieve a steady state. Conditions at the tq> of the cloud are 
now equivalent to conditions \vithin a perfectly insulated 





i.'.c .1 j;riClosure, and the radiation is ’'olacl:*'^ Tho 

i;'..''"rv nin radiation field froir, the uov/ard hemisphere incident 
c.; via-; clo.a. tep must therefore contribute to the intensity 
' ‘h.j ciroction a coiupenunt of diffusely ref iec :ed 

r .:lcn i’e. t sufficient to increase tlie intensitj'’ of radiation 
tt.^raally euittad by the cloud (in the direction ) by an 
.,..'.eunt such that the cloud appears to radiate as a biackbody-- 
fVie functional relation betv/een this diffusely scattered radiation 
and is thus the ’'mirror image" of the functional relation 

bet’veen the outgoing thermally emitted radiation at the cloud 
top and yA , and is illustrated oy the upper dashed cui’ves (IS) 
in Figures 6 and 7. • 

The lOY/er dashed curves (AS) in the figures refer to the 
same incident radiation field v/liich is now diffusely reflected? .. 
from a cloud cociposed of highly forward scattering particles 
obeying the size distribution . In this -care, on the 

ayorage, each photon is not dev-iated far. from the original 
direction of incidence even after three of four scattering pro- 
cesses , . Ducr^to the rather high finite probability of absorption 
of each phovora-particle interaction, the relative, .j^umber of .photons 
scattered back into the outward hemisphere is rather small, . 
except for those incident on the cloud top in directions near 
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grazing incidence. Again, because of the extreme forv/ard 
scattering nature of the cloud particles, even these latter 
photons must be diffusely reflected primarily xnto directions 
corresponding to very small values of ^ , Hence, the intensity 
of radiation diffusely reflected (except perhaps in directions 
corresponding jo O ) must be substantially reduced with 

respect to the intensity of diffusely reflected radiation arising 
,.i a medium of isotropically scattering centers. Increasing 
^ has the effect of raising the curves in general; 

hov/ever, we v/ould not expect a large fraction of the incident 
radiation to be diffusely reflected (due to the extreme forward 
scattering nature of the particles) until "u5q becomes almost 
unity. It sjiOuld be noted that in the limit, as the phase 
-function for single scattering becomes completely forward 
scattering [ 6> 7 ^ O \ no radiation can be 

diffusely reflected into the outward hemisphere in any direction 
C^o ^ effective emissivity of the cloud must be 

unity for all in this case. We copclude from the foregoing 
discussion that: (1) increasing has the effect of 

decreasing the overall effective emissivity, (2) the rather high 
effective emissivity at moderate to large values of yu. is due 
primarily to forward scattering, and (3) the rapid decline in 
in effective emissivity for small values oi yu. is due to departures 
from comjDlete forward scattering. 
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A-ial 02 ‘ 0 us ai’ 2 ;’ur.;ents may l>3 used to gaiu insight i:.to the 
scat'coring properties of optically thin clouds. In this case we 
■ hot:: sides of the cloud with isothermal sor.ii-inf ia:Lce 




,.ie conxriouxiOii co the *^biack racliarioa 


xhe 


V, top in the dii'ection ^ which i-esults from the diffuse 

reflpcxion of radiation originating in the upper blackbody is 
in general less than if the cloud were semi-inf initely thick, 
since some of the radiation from the upper hemisphere passes 
ccr.vpletely through the cloud in the downward direction and cannot 
be recovered as outgoing radiatio:,. An isothermal composite 


made up of an optically finite c- ud o- arlying a- semi-infinite 
perfect absorber and emitter is thus sec-i. . nave an effective 
emissivity somev/here between those of a semi-infinite cloud and 
a perfect blackbody. 

It is quite clear, the’^efore, that the sui-plus cf outgoing 
radiation at the top of the optically finite cloud models 
depicted in Figures 4a-4d and 5a-5d, compared with the relevant 
perfectly absoi-bing and emitting standards, is due to the 
temperature difference between the (warmer) ground and the cloud 
models. The difference is accentuated by the forward scattering 


nature of the cloud particles; i.e., the radiation incident on 
each cloud bottom, which is scattered, is primarily diffusely 
transmitted through the cloud in the forward direction, and 
very little is diffusely reflected back into the lower hemisphere. 
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U;;v I ‘:o Tabic 2 v/c i:ole that the effective cross 

11-^ ,') Of a unit volume mass element is smallei* 

—1 —1 —X 

T' = 5 cm"""^ and 1173.5 cm"" than at "P = 4S7.5 cm"’ 

l‘I37.3 cm -2ox the particle size distribution Of(r-j , and 
. effee :.vc albedos for single scattering, 


larger, the smaller values of — N a''^ will have the effect of 

g o ^ 

dccrv.ia:ing the normal optical thickness 7“ of a given cloud, 
and the larger values of v/ill have the effect of increasing 

the net fiov/ of outgoing radiation through the cloud. Both 
cons.tderatious v/ould lead us to expect a surplus of net outgoing 
radiation in the window region of the infrared spectrum compared 
to that which .. ^uld be calculated under the assumption of "grey*' 
clouds, where the ground is considerably warmer than the cloud. 

We should, of course, expect a deficit if the ground is cooler 
than the cloud, or if the cloud is very thick. 


7. The Net ?lux and Angular Distribution of Radiation From 
the Top of the Atmosphere 

It is evident that the outgoing Intensity Z 
the top of a cloud atmosphere is given by the relation 




- X(^i) -f 'TJkO 


( 7 . 1 ) 
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vl'.j top of tho cloud as previously considered, 'V 
is •. /. 2 troosjuission function through the ovei lying atmosphere 

is 003 direction for the v/avenumber interval of interest, 

■ / 

oof [y'‘y) is the outgoing intensity at the top of the 
overlying atmosphere due to thermal emission from that stratum. 
Loth nnd v/ere calculated v/ith the program 

furnished by F, Van Cleef v/hich v/as mentioned previously. 

Figures Sa~Sb and'9a-9b illustrate the limb darkening at 
rhe top of the atmosphere for each of the previously considered 


cloudy atmosphere models, and Table 6 gives the net outgoing 
flux 'TT P , derived from the expression (cf. Table 3) 



for each of these models. All values in the table are relative 


to the fluxes for clear sky conditions, and these fluxes have 
boon normalized to unity for each wavenumber interval. 



the 6^'c'.ucaphox’C ;:oi%ualizod to the 


I 

f— 4 
t 




% 
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io ii.;;,-.odiatcly obvious from the figur-'S that thin clouds 
cc;. jruhly cooler than the effective raaiating surface below 
. ro iv.ost effective in producing strong limb darkening. In 
. .risen with the relevant standards, limb darkening in 
.y.il-2 rases is some\vhat less at moderate values of and 

at very small values of ^ . The low cloud cases studied 

show limb darkening comparable in magnitude to those of clear 
sky conditions, while the high cloud cases show much stronger 
darkening. 

Atmospheres containing optically thick clouds still show 
limb darkening, but to a much smalle iegree. Since v/e have seen 
qualitatively that the sharp increase in darkening at sr*i all 
values of ^ is a function of the degree of forward scattering 
by individual particles, it should be possible in principle to 
obtain an idea of the effective particle sizes near the tops of 
thick cloud systems through an analysis of limb darkening curves 
obtained by airborne radiometers. Measurements near grazing 
angles would be most critical. 

Table 6 shonrs, in effect, the relative surplus (compared witu 
the relevant standai^ds, subscripted with B) of the net outgoing 
flux at the top of the atmosphere in which are imbedded clouds ox 
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i ...l/ _ optic;,! thickness, and the relative deficit of flu:>: 

'..•..tjic th^ optical thicknesses of the iubedded alouds arc infinite, 
j oeeii that the respective surpluses and deficits are most 
- V, for - 1173.5 cra"’^. In particular, the outgoing net 
r.i:: .,t the top of the atmosphere containing a semi-infinite 

'the top of v/hich is located at P =» 417 mb) is only ,927 of ’ 
t.... veutgoing net flux computed on the basis of assuming the cloud 
t-/ oj black. Conversely, the corresponding flux ratio at the 
tap of the atrAOsphere containing a cloud of normal optical thick- 
ness = 0.1 at the same level is 1.050., For a cloud of * 

normal optical thickness S, =* 0,5 this ratio becomes 1.203. Othcitif ,:j! 

"V" 

ratios derivable from the table, while not so dramatic, follow ,'ht- 

f, - 

the sa-'.ie pattern, - 



8. Couputationnl Errors ' 

it is desirable to retain as many termf^ in (3.3) as possible 
for obvious reasons (cf. Figures 3a-3d) • There will, however* be 
a poxnt beyond which N cannot be increased without destroying 
acceptable accuracy in the solutions to (4.5) and (4.20) as a 
result of rounding errors in the coxputer program. Oae critical 



area of computation centers around the simultaneous solutions of; pw 


(4.6) and (4.7) for the values of 0(s and 

^Ohi) . A check on the accuracy of tbs 

various values of vas Obtained frob tbe etxprecslba 
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5. 


2 i ' 

. i 

L 


■ A- 


- -A" 


I '}y-^' "-A' 


- (C) 


(S.l) 


-J 


:t / ; /'- <i?3 w , . 

yv - ■ ^ ^ ^ ^ ^ y ^ 


loj Cell vcilues of vl’, (ct;' ■ - ^; •••^-^'") » is (4.6) 

tho': 'ciie gi’eutest trouble will occur for jV,f | < / and = /^ ; 

this was in fact verified numerically from (S.l) , The criteriura 
f l,n.:.lly >arrived at by trial and error was that the raagnitude of 

the left-hand side of (S.l) should always be less than 10“"^ in 

Order to guarantee a comparable accuracy in (4.5) and (4.20) 
for T) :£ / . This required, under the present computational 

scheme, a limiting of the values ox N and n to those indicated 

in this paper. It should be noted that increasing the magnitude 
of T, beyond unity will correspondixigly deci-ease the 

accuracy of (4.5) and (4.20), since terms of the form Q~ ' 
occur in pairs in these equations, and the matrix solutions for 
the 2n coi atants M±oi (c<~ from the 2n boundary 

CO’-, Ij -ons [cf. eq. (4.11)] are critically dependent upon the 
difference in magnitudes of the vax'ious terms i. , the systems of 
equations to be solved. In order to be on the safe side it was 
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. cl CO rcsci’ic 
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t the study of optically finite cloua 
normal optical thicicness ^ 0,5. 

on 9, this does not unduly restrict the 


As 


we 


isi Vv ^ ^ *. k_f 1 A O S Li V< C*. y U 

..nocher important source of error centers around replacing 
th-i verious integrals encountered v/ith Gaussian sums. Since 
sin,;lo scattering contributes significantly to the diffuse 
r .di,; tion field for the values of used, the diffuse 

I’l-flection and transmission functions S and T tend to maintain 
i/.o oscillatory nature of the phase function for single scattering 
(of. Figures 3a-3d) , and are thus not amenable to reproduction 
t'-.’Cugh the use of polynomials of lew degree. Hence, the use 
v-i a quadrature formula of low degree [eq, (5,4) ] for evaluating 
lie iCitegrals in (5.3) will be source of considerable error. 

-hj maximum errors were found to be involved in calculating the 
i-,:oy:’als of the scattering function illustrated in Figure 7, 
and the integrals of the transmission functions illustrated in 
Figures 5c-5d, as indicated b 3 ^ the crosses in each figure. It 
is quite clear that the actual ’computations will not yield 
acceptable accuracy as they srand. 
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llovevev , it is extremely desirable to limit the number of 
explicit calculations of and in (5.3) 

j 

to a minimum, since each of these computations is in itself 
quite involved. The artifice of correcting the computations 
by means of physical considerations brought out in the last 
section appear to yield quite satisfactory results. The format 
for making and applying these corrections is given below. 

1. The computed values of are assumed to be 

correct as they stand.' It should be noted that the oscillatory 
nature of the phase function for single scattering is not 
maintained in the thermal component of the intensity, due to the 
diffuse nature of the radiation source. The correctness of the 
quadrature formula incorporated implicitly in the solutions of 
and thus also in (4.5) and (7.1), would appear to be 
established in view of the essentially identical results obtained 
for (>< - several test computations 

involving different values of N and n. In particular, the 
values of N tested were N = 7, 9, 11, 13, and 15, and the 
corresponding values of n were n = 4^ 5, 6, 7, and 8; only for 
N = 7 and n = 4 did the smooth curve fitted through the computed 
values differ from the smooth curves fitted through the values 


for other cases 
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2. The only additional component to the outgoing radiation 

f?Lolcl at the top of semi-inf ini je clouds is the one due to 
diffuse reflection. The relevant integral (4.13) was evaluated 
for conditions equivalent to those inside a perfectly insulated 
isothci-mal enclosure; i.e., was set equal to B(T^) 

for ail ^ ^ ~ j evaluated 

for all Cj - /j .. .^7 ) • Each computed value of 

was then corrected by the relative amount required to agree with 
the corresponding mirror image of • These 

relative corrections were then applied to the integrals computed 
which describe the intensity of radiation arising from the 
anisotropic radiation in the upward hemisphere that is diffusely 
reflected into the relevant direction , In all cases 

considered the absolute corrections required were very small, 
primarily because the diffusely reflected component its'‘'lf is 
quite small (cf. Figures 4e-4f and 5e-5f) . 

3, Corrections for the diffusely transmitted components 
through clouds of small optical thickness were computed in 
essentially the same way. In this way all the values of 

and (4.13) and (4.14) were set equal to B(T^) . 


i 
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Thc ccvaputed values of j, and i 

(5„2) were then added pnd the sum subtracted from B(T ) . The 

c 

oalculatod values of were then corrected to these 

values, and the relative corrections thus obtained applied to 
the relevant diffuse transmission integral in (5.3) describing 
the radiation diffusely transmitted through the cloud arising 
from the anisotropic radiation field of the lower hemisphere. 

The errors introduced by assuming X 5 to be exact 

should be negligible by virtue of the smallness of 
itself. The corrected points, illustrated for example in Figure 
4c-4d and 5c-5d by the open circles, were gratifyingly close 
to a smooth curve fitted through the points in all cases. 

The essential assumption required in making these correctior 
is that the form of the integrals thus corrected does not depend 
critically on the nature of the incident radiation field. Since 
only the integrals over the transmission functions would appear 
to be a cause of much concern in the cases studied, and since 

i 

the radiation field in the lower hemisphere is essentially 
isotropic in all cases, the corrections applied may confidently 
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be expected to effect a quite accurate solution. Inaccuracies 
of up to one percent may still be inherent in solutions. 

Tests of other possible sources of error in the computations 
were made wherever required, and all errors thus found were 
deteriaiaed to be negligible — on the order of 10 or less, 

9, Effects of an Arbitrary Temperature Gradient and Normal 
Optical Thickness 

The cloud models adopted in this study have been isothermal. 

Table 7 shows the effect on which is obtained 

by assuming a moist adiabatic lapse rate of 5K km“^ through a 

semi-infinite cloud. The particle densities assumed are N^ = 100 
—3 

particles cm for a cloud top temperature T = 259K, and 
N^ ^ 300 particles cm“^ for T^ = 287K. It is seen that the 
effect of a realistic temperature gradient through the cloud is 
quite small in comparison with a zero gradient, and hence 
may in general be quite safely neglected. 



Table 7, Effect of a temperature gradient ^ the cihgular </istribution of outgoing 
Radiation at the top of a 5emi- inf inite cloud 
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Tne effect of a temperature gradient on an optically 

thin cloud v/ould be expected to be much less. An idea of the 

relation among T, and calculated with the aid of equation 

(2,2) and Table 2, is illustrated in Figure 10 for the particle 

—3 

density = 100 particles cm previously cited. This figure 
is aJLso useful for obtaining an estimate of the normal optical 
thicknesses to be associated with the geometrical thicknesses 
of various clouds. 

Consider now a perfectly absorbing and emitting (nonscattering 

cloud having a temperature T and overlying an isotropically 

emitting ground radiating at twice the rate characteristic of a 

semi-infinite blackbody of temperature T . Figure 11 is a 

c 

parametric representation in terms of T, of the angular dis- 
tribution of radiation at the top I'f the cloud, computed from 
expression (5.5), where = 2B(T^) . It would appear 

from the figure that the values of t, considered in this 
study are quite representative of the extreme cases of 
interest, and ■'I’hat judicious interpolations for the scattering 
properties of clouds having other optical thicknesses are 
quite feasible. 
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10. Conclusions 

Some general summarizing statements may now be made 
regarding the infrared radiative char actex*i sties of plane- 
parallel clouds o If a cloud is much cooler than the effective 
radiating level below the cloud, it follows that; 

1. A considerable surplus flux of outgoing radiation in 
the window region of the infrared spectrum would be expected 
for thin clouds compared to what would be calculated if 
scattering were not taken into ; ccount. In general we would 

also expect clouds to be considerably poorer absorbers in the 
windov/ region than over the infrared spectrum as a ^ ::>le. This 
is due to the fact that clouds are both optically thinner and 
better scatterers in the window. We would furthermore expect 
the scattjred radiation to be predominantly diffusely trans- 
mitted through the cloud rather than diffusely reflected by 
the cloud because of the forward scattering nature of the 
individual particles making up the cloud. This in general would 
lead to a ’’surplus" of outgoing radiation in the window. 

2. A deficit of outgoing radiation would be expected In tbo 

window from optically thick clouds (i.e. £ < / ), and from 

thin clouds warmer than the effective radiating layer -below them. 
This deficit, however^ is minimized considerably by the forward 
scattering nature of individual particles. 
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3, Limb darkening in general is evident in all cases to a 
greater extent than would be anticipated by treating clouds as 
infinitely thick isothermal blackbodies. Limb darkening is 
considerable at vex'y large zenith angles, and even more greatly 
enhanced in general ii the clouds are thin and much cooler than tl 
cfxective radiating layer below them. 

4. A realistic t< mperature gradient in clouds is of aim >st n 
importance. It should be noted, however, that we have considered 
only single layers; multiple decks of clouds have not been treate 

In practice the concept of plane-parallel clouds is quite 
unrealistic except in special circumstances. Clouds are generall. 
quite ’’corrugated" at the top. These corrugations may take the 
form of mounds, billov/s, waves, and many other shapes. Furthermoi 
clouds do not have infinite horizontal extent, and because they ar 
often quite thick geometrically there is a projection factor to be 
taken into account. This projection factor is necessary to proper 
account for the apparent greater cloud cover at large zenith angl< 
than at small angles. Any general statistical analysis of the 
radiative characteristics of the Earth’s cloudy atmosphere must 
therefore be modified in accordance with these views. 

Nevertheless the preceding calculations are of considerable al 
in obtaining a qualitative picture of what one would expect to fin 
from statistical studies. In particular, if clouds are primarily 
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cumulus in type, one would expect € to be slightly less than 
unity and almost independent of zenith angle. This follows from 
considerxng cumulus clouds to be optically thick and somewhht 
hermi spherical in shape. Any zenith angle w'Ould correspond 
roughly to normal viewing incidence because of the hemispherical 
shape of the cloud. Hence, e (cf. Figures 6-7) would be rather 
high regardless of the value of . 

If thin clouds over very warm surfaces are ^resent, a surpl 
flux of radiation would still be expected in the infrarod window 
for the reasons previously cited, since very thin clouds would 
be expected to behave as though they were stratified. If the 
effective radiating level were cooler than the cloud a flux 
deficit could be expected in the window. 

Multiple cloud decks have not been treated because of the 
considerably greater complexity in the theory resulting from 
the scattering of radiation back and forth between layers. 

In reality, of course, multiple cloud systems are not uncommon. 
If the particle size distrib«.tion remains constant from layer 
to layer, and the effects of the gaseous atmosphere between 
layers is ignored, then, in the first approximation, the several 
d3cks may be regarded as one continuous cloud. Becax«se this 
"single** cloud is extended over a large vertical distance, the 
temperature becomes a strong function of optical depth. Thus 



B( C ) is no longer sensibly constant over the entire range 
of T . The effect of increasing B( "TT ) with increasing 
T would be that of causing stronger limb darkening than 
calculated previously in this paper (cf. Table 7). 

It would be very useful to extend the investigation to 
include more intervals of the infrared spectrum, as well as to 
expand the investigation to include more particle size dis- 
tributions and values of the normal optical thickn'^ss . It 

should further be n .ed that an extension to ice clouds would 
be of paramount ii forest, since absorption by ice in the infrared 
windov/ is much lower in general than absorption by water 
(KisJovskii, 1963), Consequently in general should be much 

higher for ice particler and deviations from unit emissivity 
possibly quite large for clouds composed of such particles. 

Since these particles (crystals) are in genera? not spherical, 
the Mie theory, strictly speaking, is inapplicable in describing 
single scattering; it follows that a discussion of effective 
eraissivities of ice clouds is beyond the scope of this paper. 
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APPENDIX 


Solution to Ig in the Context 

of Isotropically Scattering Particles 

The law of diffuse reflection from a semi-infinite medium 
in the context of isotropically scattering particles is given 
by (Chandrasekhar, 1960) 

where 7T^ is the incident flux of beam radiation from a point 
source which crosses (in the direction ) a unit area normal 
to the beam, and where H is defined by the nonlinear Integral 

equation 

/ 

^ f at/j' (A.2) 

By virtue of (4ol3) and (4.17) the intensity of radiation from 
all dii'ections in the upward hemisphere which is diffusely 
reflected into the direction ^ is 

/ 


(A. 3) 
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or, from (A.l), 






jMo 




(A. 4) 


Now let 


(A. 5) 


for all , where is the Planck function over the 

interval Av characteristic of the cloud temperature T^, This 
artifice renders conditions a' e cloud top equivalent to 
conditions within a perfectly insulated isothermal cavity, and 
(A. 4) becomes 




(A.6) 


Upon using (A. 2) 
with 


and replacing 


in the integrand of (A.6) 
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Again, upon multiplying (A. 2) by and integrating over 

the range of ^ , we obtain 




> 


(A. 9) 


Or, upon interchanging ^ and and adding, 




~ +- 



(A. 10) 


Solving the quadratic in (A, 10) for the intergral gives us 

/ ‘ ^ 

^^0 - ) ± , (A. 11) 

O 

Since /^) converges uniformly to unity as uniformly approaches 
zero [cf, eq. (A. 2)], we must choose the minus sign in (A. 11) in 
order to preserve the identity for all , 


By virtue of (A. 11) we obtain for (A. 8) the result 



(A. 12) 
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V/e have solved the problem in the context of conditions 
within an isothermal cavity, where the cloud top is one of the 
perfectly insulated bounding surfaces. It follows from 
Kirchhoff's law that the thermally emitted component of radiation 
from the cloud is given by 

m 

or, from (A. 12), 

BjJti) H( 6 <) y (A. 14; 

and the desired result is established. 



REFERENCES 


Chandrasekhar, S,, 1946 : On the Radiative Equilibrium of 

a Stellar Atmosphere. XII, Astrophys, J «, 1^, 191-202, 

, 1960 : Radiative Transfer . New York, 

Dover Publications, Inc,, 393 pp. 

King-, Jean I. F., and E, T. Florence, 1964 ; Moment Method for 
Solution of the Schwarzschild-Milne Integral Equation. 
Astrophys. J, , 339 , 397-403. 

Kislovskii, L, D.^ 1959 : Optical Characteristics of Water and 
Ice in the Infrared and Radiowave Regions of the Spectrum. 

Opt, i Spec ., 7^, 201-206. 

Neiburger, M., 1949 : Reflection, Absorption, and Transmission 
of Insolation by Stratus Cloud. J. Met ., 6, 98-104. 

Penndorf, R. B. , 1963 : Research on Aerosol Scattering in the 
Infrared, Final Report AFCRL-63-668. Published as 
Research and Advanced Development Division AVCO Corporation 
Technical Report, RAD-TR-63-26, 231 pp. Available at cost 
through the U, S, Department of Commerce, Office of 
Technical Services, Washington , D. C. 20235. 

Saiauelson, R. £., 1965 : Radiative Transfer in a Cloudy Atmosphere 
NASA Technical Report, R-215. Available at cost through the 
U. S. Department of Commerce, Office of Technical ’ Services, 
Washington , D. C. 20235. 



REFERENCES (Continued) 


Wark, D. Q,, J. Alishouse, and G. Yamamoto, 1964 : Variation 
of the Infrared Spectral Radiance near the Limb of the 
Earth, Appl, Opt ., 221-227. 

, G. Yamamoto, and J, Lienesch, 1962 : 

Infrared Flux and Surface Temperature Determinations from 
TIROS Radiometer Measurements. U. S. Weather Bureau, 
Meteorological Satellite Laboratory Report No. 19, 84 pp. 
Available at cost through the U. S. Department of 
Commerce, Office of Technical Services, Washington , 

D. C. 20235, N-63 - 11881. 



FIGURE LEGENDS 


Figui'e 1. Complex index of refraction of water (^= yt ^ k ) 
as a function of the wavenumber 1 / in inverse 
centimeters. 

Figure 2. Particle size distributions (modified from Nc»uurger) 
representative of California stratus and the Mie 
efficiency factors Q 2 and as a function of 1 / and 
particle diameter d. The computed points for Qg and 
Qa are indicated respectively by circles and triangles 
Figure 3. Phase functions for a single scattering /»(cos ^ ) in 
arbitrary units as a function of the scattering angle 
Figure 4. Contribution to the total intensity (C) at the cloud 
top from the components respectively diffusely 
reflected (S) , themally emitted (E) ^ directly 
transmitted (D) , and diffusely transmitted (T) for 
the wavenumber interval centered about V == 889.5 cm"”^ 
All Intensities, including the standard of compariscMQ 
(B) , are on the same (but arbitrary) scale. 

Figure 5. Same as Figure 4 for the wavenumber interval centered 
about y = 1173.5 cw“^. 



FIGURE LEGENDS (Continued) 


Figure 


Figure 

Figure 


Figure 


Figure 


Figure 


6. Effective emissivities (solid curves) *and diffuse 
reflectances (dotted curves) of semi-infinite clouds 
for a>o = .403. The curves labeled (AS) and (IS) 
refer to anisotropic and isotropic scattering 
respectively. All curves are normalized tc unity. 

7. Same as Figure 6 for ui « .768. 

8. Limb darkening at the top of the atmosphere in the 
wavenumber intei'val centered about y = 889.5 cm”^. 

The normal 'optical thickness of each cloud model 
imbedded in the atmosphere is indicated by the 
relevant value of T, , and each cloud is located 

at a height characteristic of the indicated pressure P. 

9. Same as Figure 8 for the wavenumber interval centered 
about V = 1173.5 cm*"^. 

10. Relation among the temperature difference aT^ 

gecmietrical depth AZ, and normal optical depth T for 

—3 

a particle density 100 particles cm and a 

cloud top temperature = 259K. 

11. Limb darkening for isothermal nonscattering clouds of 
temperature as a function of the parameter T, • The 
ground is assumed to be a blackbody of temperature 
radiating at a rate B(T^) » 2B(T^) • 

















o 


4.0 

















o 



d (MICRONS) 



6 (degrees) 





















2 


COMPUTED CO 
CORRECTED (T) 



















EMISSIVITY 





1 ( 0 . 1 ) 


SCATTCftlNGMOOCLS 
NON-SCA7TEIUNC STANDARDS 


J I L 


.0 


1.0 




7 











